Probing single-photon ionization on the attosecond time scale 
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We study photoionization of argon atoms excited by attosecond pulses using an interferometric 
measurement technique. We measure the difference in time delays between electrons emitted from 
the 3s^ and from the shell, at different excitation energies ranging from 32 to 42 eV. The 
determination of single photoemission time delays requires to take into account the measurement 
process, involving the interaction with a probing infrared field. This contribution can be estimated 
using an universal formula and is found to account for a substantial fraction of the measured delay. 



The interaction of light with matter is an essential pro- 
cess in nature and its paradigm, the photoelectric effect, 
has been studied during decades using synchrotron ra- 
diation [1 . The development of ultrashort light pulses 
in the attosecond range allows scientists to tackle tem- 
poral aspects of electron transitions in atoms, molecules 
and more complex systems. Cavalier i et al. [2] investi- 
gated photoemission from the valence and the conduction 
band in tungsten crystals using single attosecond pulses 
and an infrared probing field through the streaking tech- 
nique [S', T. Recently, Schultze et al. [5 implemented 
the same technique to study photoemission from the 2s^ 
and 2p^ shells in neon at a pulse energy of 100 eV. They 
measured a difference in photoemission time delays equal 
to 21 as, a value which is significantly larger than the ex- 
pected theoretical value, as further discussed in a series 
of theoretical articles |6tt9j. 

In this letter, we examine photoemission of electrons from 
the and 3p^ shells in argon. Our method uses a fre- 
quency comb of high-order harmonics with photon en- 
ergies varying from 32 to 42 eV for the photoionization 
and a weak infrared laser field for probing the outgoing 
electrons in time. Our method is complementary to the 
streaking method used in [5 . It is based upon interfer- 
ometry and it allows us to explore the threshold region 
for the 35^ shell, where one expects large variation in 
photoemission times. We measure a delay between the 
ionization from the 3^^ and 3p^ shells which varies with 
photon energy. We investigate the influence of the in- 
teraction with the weak infrared field, which is needed 
to get the interferometric measurement and therefore the 
temporal information. Probing the outgoing wave packet 
with even a weak infrared field affects the electron mo- 
tion and therefore the measured delay. Fortunately, this 
effect can be analytically calculated and takes a universal 
form, which allows us to disentangle the different effects 
and gives us access to the single-photon ionization time, 
also called Wigner time [lOl [11] . 

The basic principle of our experiment is shown in Fig. [l] 
We ionize argon using a comb of high-order harmonics. 



For a central frequency of the harmonic comb above the 
binding energy of the 35 orbital we simultaneously create 
two independent electron wave packets, one originating 
from the 35^ and one from the 3p^ subshell. The presence 
of a fraction of the fundamental laser field with frequency 
uj induces the formation of sideband peaks due to two- 
photon transitions including absorption or emission of an 
infrared photon [121 US]- Two different and interfering 
quantum paths involving consecutive harmonics lead to 
the same sideband (see Fig.[T]). When changing the delay 
r between the harmonic comb and the laser field, the 
sideband signal from a given subshell is modulated as \\A\ 



S{t) =a + /3cos[2cj(r + rA + ri)], 



(1) 



where a, j3 are two constants independent of r. The term 
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FIG. 1. (color online). Principle of the measurement. Two 
electron wave packets originating from different subshells are 
simultaneously created using the same comb of high order 
harmonics. The outgoing electron wave packets are further 
probed with a weak infrared field. For simplicity only two 
harmonics are indicated. Also shown is the experimental har- 
monic spectrum used. 
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Ta is proportional to the difference in phase between 
consecutive harmonics and describes the group delay 
of the attosecond pulses, while Ti represents the atomic 
delay due to the two-photon ionization process [T5HT7] . 
As we will show below Ti can be connected to the 
Wigner time delay Tw for the single-photon ionization. 
The knowledge of Ta as well as of the absolute value of 
the delay r would enable us to determine Ti directly. 
However, these are difficult to obtain experimentally. 
The simultaneous measurement of the two electron wave 
packets allows us to cancel the influence of the attosec- 
ond group delay Ta and to determine ri(3s) - ri(3p) 
at the same photon energy, i.e., at kinetic energies 
separated by the difference in binding energy between 
the two orbit als (13.5 eV). 

Our experiments were performed with a 800-nm 30- 
fs Titanium- Sapphire laser system [18^. High-order 
harmonics were generated in a pulsed Ar gas cell and 
spatially filtered using a small aperture [19]. We used a 
0.2 /im thick chromium thin film to select a 10 eV-broad 
spectral window corresponding to harmonic 21 to 27 
at 38 eV central energy (see also Fig. [T]). This filter 
was chosen to separate the wave packets emitted from 
the 35 and 3p subshells in energy. The comb of about 
four phase-locked harmonics, corresponding to a train 
of attosecond pulses with a duration of 450 as, was 
focused by a toroidal mirror into the sensitive region 
of a magnetic bottle electron spectrometer containing a 
diffusive Ar gas jet. Part of the laser field was extracted 
prior to the high-order harmonic generation and recom- 
bined downwards collinearly with the harmonics with a 
variable time delay r. The precision of our measurement 
does not depend on the duration of the attosecond 
XUV pulses but on the interferometric stability of our 
experiment. 

Fig. |2|a) and 2(b) present electron spectra as a func- 
tion of the delay r between the XUV and the infrared 
pulses. The low-energy spectrum in Fig. [2|a) shows elec- 
tron peaks at energies corresponding to single photon ion- 
ization from the 35 shell by the harmonics and additional 
sideband peaks due to two-photon transitions. The high- 
energy part of the spectrum shown in Fig. |2jb) presents 
the corresponding photoelectron spectra for 3p ioniza- 
tion. Although simultaneously recorded the results are 
presented separately due to the unequal signal strength 
caused by the difference in cross section and detector sen- 
sitivity (note the different color scales). For both chan- 
nels the sideband signal oscillates. Fig. [2|c) presents the 
delays obtained by Fourier transform of the sideband sig- 
nal along the delay axis for the scan shown in (a) and 
(b), corrected for the influence of the Cr filter, which is 
positively dispersive in this region [20 . The variation in 
delay reflects mainly the positive chirp of the attosecond 
pulses. The main experimental result of the present work 
is the significant offset between the delays measured for 
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FIG. 2. (color online). Energy spectra as a function of delay 
from electrons liberated from the 3s orbital (a) and the 3p 
orbital (b), respectively, (c) Retrieved delays corrected for 
the Cr group delay. Also shown are the 3p-delays shifted 
down in energy for comparison with the 3s-delays (dashed 
line) . 



the two wave packets. To emphasize this result, we show 
as a dashed line the 3p-delays shifted down in energy by 
13.5 eV. Taking the difference between the measured de- 
lays at the same excitation energy and averaging over five 
independent measurements, we determine a difference in 
delays ri{3s) -Ti{3p) equal to -40±10 as for sideband 22, 
-llOilO as for sideband 24 and -80±30 as for sideband 26. 
To understand the meaning of these time delays, we need 
to establish the connection between single-photon ion- 
ization and the two-photon ionization process used in 
the measurement. The phase of the transition matrix 
element describing a single ionization process towards a 
final state with angular momentum £ is the scattering 
phase r]£^ i.e. the phase accumulated by the photoelec- 
tron when escaping from the atom. Its energy derivative, 
Tw = hdr]£{e) jde represents the "photoionization time de- 
lay" also called Wigner time delay [TOllTl]. Clearly, both 
and Tw depend on the details of the atomic potential 
and their computation remains a challenge for theory. 
Using second-order perturbation theory, the transition 
matrix element for two-photon ionization involving ab- 
sorption of a harmonic photon cJh and a laser photon uj 
from an initial state Lpi to a continuum state Lp^ with 
asymptotic momentum k can be written as 



-iEjEii lim 

£^0 + 



((Pj^\e-f\(pn)((Pn\e-r\(pi) 
ei + cJh - en + is 



. (2) 



Atomic units are used throughout. The complex ampli- 
tudes of the laser and harmonic fields are denoted E^^ and 
Eii and e is their common polarization vector. The en- 
ergies of the initial and intermediate states are denoted 
Ci and e^,, respectively. The integral-sum is performed 
over all possible intermediate states cpn- The index a in- 
dicates that we first discuss a two-photon process with 
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absorption of the laser photon. 

We consider the channels s ^ p ^ £ with £ = s^d. Us- 
ing spherical coordinates, separating radial and angular 
parts of the wave functions, and expanding the final wave 
function into partial waves, the transition matrix element 
becomes 

Mf )(fc) = -iE,E, ^ CeoYeoik) e'^'^'^^T^^^k) , (3) 

£=0,2 

where Y^o is a spherical harmonic, C^o the correspond- 
ing angular coefficient, and rji the scattering phase of the 
final state. The radial two-photon transition matrix ele- 
ment T^\k) can be expressed as p!3l [T5] 



T^'\k)=i 

Jn 



{RM\r\Rni){Rni\r\RiQ) 



{RM\r\pk^i). (4) 



In the right member of equation Q we introduce the 
perturbed wave function p^^i with the wave number ka 
such that kll2 = + = k'^ 1^ - uj (see Fig. [l])[2r. To 

(2) 

get an estimate of the phase of , we consider the 
asymptotic behavior of the wave functions involved in 
equation Q. The perturbed wave function pk^i is an 
outgoing wave [22l [23] 

1. /X i\kar-^T^+^\n{2kar)+r)i{ka)\ 

\im pk^i{r) oceV ^ -I, (5) 

r— >cxD 

while Rki is real with an asymptotic behavior: 

lim Rkiir) oc sin [kr-ln+l\n(2kr)+m(k)]. (6) 

The factor £7t/2 arises from the centrifugal potential, 
while ln{2kr)/k is a correction due to the long range 
Coulomb potential. Using Eqs. (Isl) - (IgI) we find an ap- 
proximate expression for Ma (k) 



(I) 



\ka-k) 



{2k)k 

i^kaW^ 



r(2 + zz), (7) 



(11) 



where z = 1/ka - l/k and T{z) is the complex Gamma 
function. The first phase term (I) is the scattering phase 
of the intermediate state and identical to the phase of 
the corresponding one-photon ionization. The phase of 
term (II) can be assigned to the laser-driven transition 
connecting the two continuum states in the presence of 
the long-range Coulomb potential, (p^^. It is independent 
of the short range behavior of the atomic potential and 
therefore universal. Corrections to this approximation 
due to the core are expected to become important only 
at energies close to threshold. 

The phase of the two-photon matrix element M^'^^ for the 
second pathway, i.e. absorption of an harmonic photon 
cJh followed by emission of an infrared photon uj via an 
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FIG. 3. (color online). Computed delays associated with the 
following ionization channels: (a) 3p es, (b) 3p ed, (c) 
3s ep in Ar, and (d) Is ep in H. The red lines are the 
one-photon Wigner time delays. The blue lines represent the 
estimated delays induced by the measurement, Tcc- The sum 
of the two delays is shown as a black line. The dashed line in 
(d) is the result of an exact calculation in H. 



intermediate state with wave number kl/2 = /c^/2 + co' (see 
Fig. [T]), can be derived in a similar manner. The total 
interference signal is obtained by angular integration of 



Mi^^ + Me^^^l^. It can be written as Eq. with 



r(2)|2 



2u 



2uj 



(8) 



This result gives an intuitive understanding of the ion- 
ization time Ti{Zs). It can be expressed as the sum of the 
Wigner time delay Tw for one-photon ionization 3s ep 
and an additional continuum-continuum delay Tcc inher- 
ent to the measuring process. This analytical derivation 
can be easily generalized to other ionization channels. 
Fig. |3] shows the delays involved in the three ionization 
channels 3p es (a), 3p ed (b) and 3s ep (c) in Ar as 
a function of kinetic energy. The Wigner time delay Tw 
(red) is obtained by taking the derivative of the scatter- 
ing phase taken from [24]. For comparison, we also show 
in (d) the delays for the pathway Is ^ in hydrogen in 
the same energy region. The continuum-continuum de- 
lay Tcc (blue) is calculated for a 800 nm laser wavelength 
and identical for all the channels and atoms. The black 
line indicates Ti as the sum of the two contributions. The 
Wigner time delay variation can be nicely and intuitively 
interpreted. Low-energy electrons take longer time to es- 
cape from a given subshell than high-energy electrons. 
Furthermore, electrons escaping to a channel with higher 
angular momentum take longer time than those escaping 
towards a channel with low angular momentum because 
of the centrifugal barrier. The continuum-continuum de- 
lay has the opposite behavior and leads to an apparent 
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FIG. 4. Comparison between the measured delay differences 
for ionization of Ar from the 3s and 3p shells (diamonds) with 
calculations performed according to the approximate theory 
developed in this work (black line). Also shown is the delay 
expected for one-photon ionization (red line) and the laser 
driven continuum-continuum transition (blue line). 



sented here will stimulate numerous experiments, aiming 
at measuring photoemission delays in relative or possibly 
absolute values, in a variety of systems. The recent de- 
velopment of high-order harmonic generation using mid- 
infrared drivers [26, 27 should lead both to a higher spec- 
tral resolution for this type of measurements as well as 
to a broader investigation bandwidth. 
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supported by the European Research Council (ALMA), 
the Marie Curie Intra-European Fellowship ATTOCO, 
the Knut and Alice Wallenberg Foundation, the Joint Re- 
search Programme ALADIN of Laserlab-Europe II and 
the Swedish Research Council. We acknowledge financial 
support from the ANR ATTO-WAVE. P.J. acknowledges 
support from the Swedish Foundation for Strategic Re- 
search. 



quicker escape for the low-energy electrons. Finally, we 
also indicate in Fig. [sj^d) results from exact calculations 
in H (dashed line). The comparison between the black 
and dashed lines gives an estimation of the error made 
in considering only the asymptotic behaviors of the per- 
turbed and final wave functions. 

For the energy range considered in the present work the 
asymmetry parameter remains close to two [25 , which 
indicates that the ionization channel 3p ed dominates 
over 3p ^ es. Neglecting the 3p es channel, we calcu- 
late ri(3s) - Ti{3p) at the same excitation energy, using 
the approximation presented in Eq. (|8| . Fig. [4] presents 
the approximated delays (black line), together with the 
experimental results (o). We also show the Wigner time 
delay (red line) and the influence of the laser-driven 
continuum-continuum transition (blue line). The exper- 
imental results at the two highest energies agree well 
with the results of our calculation. The lowest energy 
point, however, lies outside our estimated uncertainty. In 
this region the core plays a more important role for the 
continuum-continuum transition, as well as the Wigner 
time delays may differ from those calculated in [24 . Us- 
ing our estimated continuum-continuum delays, we can 
tentatively deduce the difference in photoemission delays 
to be equal to 140 as at 34 eV and -20 as at 37 and 40 eV. 
In conclusion, we have performed experimental measure- 
ments of photoemission from the 3s^ and 3p^ shells in 
Ar, using a weak laser field to probe the created electron 
wave packets by interferometry. We identify two contri- 
butions to the measured delays: the Wigner time delay, 
which represents the group delay of the electron wave 
packet emitted in single photon ionization, and a delay 
inherent to the measurement process. It can be simply 
estimated using an universal formula which only depends 
on the laser frequency. We believe that the work pre- 
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